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Abstract 

We present a toy model study of the high temperature deconfining transition 
in Yang-Mills theory as a breakdown of the confinement condition proposed by 
Kugo and Ojima. Our toy model is a kind of topological field theory obtained 
from the Yang-Mills theory by taking the limit of vanishing gauge coupling con- 
stant (?ym — ► 0, and therefore the gauge field is constrained to the pure-gauge 
configuration = g^d^g. At zero temperature this model has been known to 
satisfy the confinement condition of Kugo and Ojima which requires the absence 
of the massless Nambu-Goldstone-like mode coupled to the BRST-exact color cur- 
rent. In the finite temperature case based on the real-time formalism, our model 
in 3+1 dimensions is reduced, by the Parisi-Sourlas mechanism, to the "sum" of 
chiral models in 1+1 dimensions with various boundary conditions of the group 
element g(t,x) at the ends of the time contour. We analyze the effective potential 
of the SU(2) model and find that the deconfining transition in fact occurs due to 
the contribution of the sectors with non-periodic boundary conditions. 
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1 Introduction 



Deconfining transition in high temperature QCD has been a subject of active researches for 
both theoretical and phenomenological interests. This transition may be stated as follows: 
below a certain critical temperature f3~ l we cannot observe color non-singlet excitations, while 
above colored excitations such as quarks are allowed as physical ones. Widely adopted as a 
criterion of (de) confinement at finite temperature is the expectation value (P)g of the Polyakov 
loop (thermal Wilson loop) operator P(x) = trPexp (^J^drA^r, x)^j which measures the free- 
energy of a quark put in the system as an external source 0. However, the Polyakov loop 
cannot be used as a criterion of confinement in a system with dynamical (quantized) quark 
fields. Moreover, (P)p tells nothing about the (de) confinement of other color carrying fields, 
for example, the gluon. In this paper we shall study the deconfining transition on the basis 
of another confinement criterion proposed by Kugo and Ojima (KO) ||. In contrast with the 
Polyakov loop, the confinement mechanism of KO treats directly the confinement of quantized 
colored fields (particles) of any kind. 

The KO confinement criterion is based on the BRST quantized Yang-Mills theory || 
described by the lagrangianF] 

^YM = r~ 2 — tr-^/L + ^matter — i5 B G , (1.1) 

z 9ym 

where 5 B is the BRST transformation defined as usual by 

5 B A^ = D„c = <9 M c+ [A^c], 5 B c=-^{c,c}, 5 B c = iB, 5 B B = 0. (1.2) 

The last term of eq. ( |1 . 1| ) gives the gauge-fixing and the corresponding ghost terms. The key 
quantity in the KO confinement mechanism at zero temperature is the BRST-exact conserved 
color current iV M , 

N lt = -iS B K ll = {Q B ,K li }, (1.3) 

where Q B is the BRST charge, and is obtained from £ym of eq. ( |1.1[) by making a local 
gauge transformation S e , = D^e and 5 e (f> = [4>, e] for (f> = c,c, B, as 

5Xym = -tS e 5 B G = -i8 B K£ ■ d^. (1.4) 

* We restrict the gauge group to SU(N). The field variables <f> — A^, c,c,B are Lie algebra valued and 

are expressed as <ft — X^T 1 ( t ) °'^ a * n terms of Hermitian fields (j) a and the (anti-Hermitian) basis t a with the 
normalization tr(t a t b ) = ~{l/2)6 ab . 
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Note that S e £.YM has its contribution only from the gauge-term 5-qG, and the last expression is 
because we have assumed that the gauge-fixing function G preserves the global color rotation 
symmetry. In the Feynman-type gauge with G = tr [c (d^A^ — aB)], we have = D^c. The 
ordinary Noether color current J® = qjnt a q + . . . containing the matter fields is related to 
by J M = N^ — {1/ gYyi)d u F VIX using the equation of motion. Therefore both iV M and J M generate 
the same global color rotation. 

The KO confinement condition requires that the BRST-exact color current ( |1.3|) con- 
tains no (Nambu-Goldstone-like) massless one-particle component. If this condition is satis- 
fied, then the integration / d 3 xN has a well-defined meaning and hence the color charge Q a 
can be written in a BRST-exact form, Q a = {Qb, J^xKq}, which implies color confinement 
in the sense that any color non-singlet asymptotic state is necessarily BRST unphysical and 
unobservable.Q 

In order to extend the KO confinement mechanism to the finite temperature case, we shall 
recapitulate the elements of finite temperature gauge theory. First, in statistical mechanics 
of gauge theories in the BRST formalism, the statistical average must be taken only over 
physical states. Let Pphys be the projection operator to the subspace of physical states 7i p h ys = 
Ker Qb/IhiQb- Then we have the following useful identity || for the thermal expectation value 
(0)p of a BRST invariant operator O satisfying [Qb, O] = 0: 

(0) p = Tr (V phys e-P H 0) /Z((3) = Tr (e~^ N ^o) /Z(fi), (1.5) 

where iV g h is the ghost number operator, and Tr means the trace operation over all (physical 
as well as unphysical) states: 

TrO = J2(k\0\l) V rk Ui = (k\l), T,VkiV^ = S km ). (1.6) 

k,l V I ' 

The partition function Z(J3) itself is also given by Z(J3) = Tr (Pphyse-^^) = Tr L-fiB+iKN^y 
Eq. ( |1.5|) is a consequence of the formula 1 = V p hy S + {QB, 3 R} where 1 is the identity operator. 
We adopt the last expression of eq. (|1.5| ) with statistical weight e -P H + tnN gh as the definition 
of (G)a even when O is not a BRST invariant quantity. 

A framework of finite temperature field theory which is suitable for discussing the KO 
mechanism is the real-time formalism where we can treat fields with ordinary time variable 
t (— oo < t < oo). There are two formulations of the real-time formalism. One is the path- 
integral formalism pi, and the other is the operator formalism called thermo field dynamics 

t Detailed explanation of the KO confinement mechanism at zero temperature is found in the original paper 
(U (see also ref. j| for a brief explanation) . 
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(TFD) |7|, ||. The path- integral formalism is convenient for concrete calculations, while TFD 
is necessary to generalize the KO confinement mechanism, which has been formulated in the 
BRST operator formalism, to the finite temperature case. In this paper we shall use both 
of the two formalisms regarding them as equivalent. We use the same symbol for both an 
operator and the corresponding path-integration variable. 
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Figure 1: The time contour C. 

First the path-integral formalism applied to a gauge theory is as follows. We consider 
thermal Green's functions generated by 



Z[j] = Tr { e -/3*+-iV ghTc | exp i Jj^ | ? (L7) 



where we have used the abbreviation J c A = J c dr J d 3 xA(r, x) with the contour C in the 
complex time-plane depicted in fig. [I], and Tc denotes the ordering along the contour C. <fi 
and j represents a generic field operator in the system and the corresponding source. The 
generating functional Z[j] (|L7l ) has a path-integral expression 0: 



Z[j]= J P0exp|ijT (£ym(0)+J0) 



periodic 

where the functional integration T>cf) = T>A^ VcVcVB should be done with periodic boundary 
condition 0(— T— if3) = 0(— T) for all the fields <fi = A^, c, c, B. In particular, the effect of the 
factor e lnNsh in eq. (|1.7| ) is to turn the boundary condition of the fermionic fields c and c to the 
periodic one. Since the contour C of fig. [I] contains the horizontal (real-time) segments C\ and 
C*2 of infinite length (we take the limit T — > oo), we can consider thermal Green's functions 
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with ordinary real time arguments in contrast to the imaginary-time formalism where the 
contour C is simply a straight vertical line [0, — 

Thermo field dynamics (TFD) |7|, [8| is the operator formalism which reproduces the thermal 
average of eq. ( |1.5| ) as an expectation value with respect to "temperature dependent vacuum" 
10(0)): 

{0) p = W)\OW))- (1-9) 

Application of TFD to gauge theories is given in ref. ||. Corresponding to two horizontal 
segments G\ and C 2 of the contour C in fig. [I], we have to double the fields in TFD as 
compared to the theory at zero temperature. The fields and states corresponding to C\ is 
denoted as before by (fi(t,x) and \k), and those corresponding to C2 is denoted with tilde; 
4>(t,x) and \k). Then |0(/3)) is given as 

One can easily see that eq. (|1.9Q holds for this |0 (/?)). The lagrangian £ym describing TFD is 

^YM = ^YM _ ^YM, (1-11) 

where £ym is the lagrangian for the tilde fields 4>. The physical states in TFD are specified 
by Qb Iphys) = using the BRST charge of the whole system Q B = Qb — Qb- See refs. |9| for 
precise definitions. 

Now we are ready to generalize the KO (de) confinement criterion to the finite temperature 
Yang-Mills theory. In the finite temperature case based on TFD, the relevant BRST-exact 
color current has contributions from both the non-tilde and tilde fields (cf. eq. (|1 . 1 1|) ) : 

A^ = iV M -AT = {Q B ,Z M }, (1.12) 

where K ^ = — K^. Note that the vacuum |0(/9)) (|1.10| ) is not invariant under a separate 
color rotation on non-tilde or tilde fields generated by or iV M . Using this N^, the KO con- 
finement mechanism is generalized to the finite temperature Yang-Mills system as follows. // 



the BRST-exact color current iV M contains no massless one-particle component, it im- 

plies color confinement: the system contains no colored excitations (quasi-particles) as physical 
ones.0 

In order to study the deconfining transition, we have to first prepare a situation where 
the confinement is realized at zero temperature {(3 = 00). Although the KO confinement 
* We are assuming that the asymptotic field like analysis applies also to TFD. 
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criterion has not been shown to hold in real QCD at j3 — oo, we have a toy model of four- 
dimensional Yang-Mills system where the confinement condition of KO is known to be satisfied 
at f3 — oo. This toy model is obtained from the real Yang- Mills system (|1 . 1|) by taking the 
limit of vanishing gauge coupling constant; (?ym — ► 0. In this limit, the field strength term, 
(1/ gY M )trF^ u , of eq. forces the gauge field constrained to the pure-gauge configuration, 
A^x) = g jf (x)d fJ ,g(x), and hence the system is reduced to something like a topological field 
theory described by the BRST-exact lagrangian —i8-BG[Ap = g^d^g] [H| . We call this toy 
model the pure-gauge model. 

Although the pure-gauge model contains no physical degrees of freedom, the KO confine- 
ment criterion is still a non-trivial dynamical problem. In naive perturbation theory of both 
the real Yang-Mills theory and the pure-gauge model at zero temperature, we have a massless 
pole coupled to N^. This is revealed by the Green's function, 

dAp -e^(YN a (x)A b M) p=00 ~ 5 ab ^. (1.13) 



(2vr) 4 " v ' " v y/p -°° p 

However, in the case of pure-gauge model with a special gauge called OSp(4/2) symmetric 
gauge, we can show the massless pole in (|1.13|) is in fact missing. This is because the pure-gauge 



model in 3+1 dimensions is shown to be "equivalent" by the Parisi-Sourlas mechanism |L2 
to the chiral model in 1+1 dimensions where the vacuum is realized in disordered phase with 
a mass gap ||13|| . It is expected that, in the real Yang-Mills theory also, the KO confinement 
condition (at j3 = oo) holds due to large gauge field fluctuation in the direction of local gauge 
transformation! 

Having completed the preparation, let us turn to the explanation of our analysis of finite 
temperature case carried out in this paper. Recalling that the pure-gauge model is obtained 
as the gyu —> limit of Yang- Mills theory ( |1 . 1| ) , we consider the two-point function 



in the limit (?ym ~~ * 0. (In eq. ( p.. 14 ), T denotes the time-ordering in TFD and the T^-ordering 



in the path- integral formalism.) We know that (|1.14j) is free from massless poles at zero 
temperature (5 — oo, and study whether a massless pole is generated at high temperature 
(3 < (3 C with some critical f3 c . For this purpose we observe that in the path-integral formalism 
Q1.8Q the A^-integration is reduced in the c/ym — ► limit to the ^-integration (A M = g^d^g) 



§ In fact, the color confinement by the KO mechanism is interpretable as a consequence of the restoration 
of local gauge symmetry with transformation parameter e(x) ~ a M x M (a^: const.) and hence S^A^ ~ a M 



5 



and, in addition, that, although = g^d^g has to satisfy the periodic boundary condition 
A^(—T—i(3,x) = A^—TjX), the group element g(t,x) need not be strictly periodic. In 
fact, the boundary conditions using the constant SU(N) left-transformations, g(—T—i(3, x) = 
h ■ g(—T, x), are allowed ones which respect the periodicity of g^d^g. Adopting the OSp(4/2) 
symmetric gauge of ref. JXX|] we find that the two-point function ( |1.14| ) in 3+1 dimensions is 



equal to the average of the Green's function JcPp/(2ir) 2 e ip - x (TAl(x)Al(0))p over the boundary 
conditions of g(t,x) in the 1+1 dimensional chiral model. This implies that color confinement 
by the KO mechanism breaks down if the corresponding chiral models in 1+1 dimensions are 
realized in the ordered Nambu-Goldstone phase for a finite range of the boundary conditions. 

Since we do not have a systematic non-perturbative methods to analyze finite temperature 
chiral model in the real-time formalism with unusual boundary conditions, we shall carry out 
the analysis of the effective potential of the 0(4) non-linear cr-model (which is equivalent 
to the 577(2) chiral model) obtained by the large- N like method. Our analysis shows that 
there is indeed a desired deconfining transition at 1//3 ~ m (m: mass gap at (3 = oo). The 
transition occurs because the infrared singularity which realized the disordered phase at zero 
temperature is weakened by the non-periodic boundary conditions. 

The rest of this paper is organized as follows. In Sec. 2, we introduce the finite temperature 
pure-gauge model in the real-time formalism as the (?ym — > limit of the Yang-Mills theory, 
and explain the Par isi- Sour las reduction to the two-dimensional chiral model. In Sec. 3, we 
carry out the analysis of the 0(4) non-linear cr-model with generalized boundary conditions. 
The final section (Sec. 4) is devoted to summary and discussion. 

2 The pure-gauge model 

As stated in the Introduction we shall consider the g\M —* limit in finite temperature Yang- 
Mills system in the real-time formalism described by the path-integral (O). In this limit the 
gauge field A^ is restricted to the pure-gauge configuration, 



A^x) = g\x)d^{x) with g(x)eSU{N), (2A] 



due to the (1/gy^trF 2 ^ term in C YM (|1.1[) , and the system is reduced to the pure-gauge model 
(PGM) with dynamical variables (g, c, c, B) just as in the zero-temperature case [JTOj, p| . What 
is particular to the finite temperature case is the boundary condition of g(t,x). Although 
A^ = g^d^g has to satisfy the periodic boundary condition, 

A li (-T-i/3,x) = A li (-T,x), (2.2) 
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in eq. (|1.8|) , the group element g(t,x) need not be strictly periodic. Boundary conditions 
related by a constant SU(N) left-transformation, 



B h : g(—T— i/3, x) = h-g(-T, x) (h £ SU(N)) , (2.3) 

are allowed ones which respect the periodicity of = g^d^g. 

Therefore it is natural to assume that in the limit gywi — * the A^-integration with 
periodic boundary condition is reduced to the g-integration using all the boundary conditions 
Bh] namely we have to integrate over the boundary condition parameter 

fvAp - JdhJ B Vg, (2.4) 

periodic 

where J dh denotes the Haar measure of SU(N). Then the expectation value of an operator 

is reduced in the #ym — > limit to 

(0) p -> Jdh VgVcDcDB O exp (i J^ £ PGM ) j JdhZ h , (2.5) 

where £pgm is the lagrangian of PGM, 

£pgm = — ^'&G\ Aii=g ^ diig , (2.6) 

and Zh is the partition function of the PGM with the boundary condition Bh- The BRST 
transformation on g is given by S^g = gc. The partition function Z% is shown to be equal to 

1 for any h by reversing the manipulation of eq. (|1.5|) : 



Z h = Tr (u h e^ H+i ^) = Tr (u h V phys e^ H ) = 1, (2.7) 



where Uh is the operator of the left-transformation by h, U h ~ 1 gUh = h-g. In eq. ( |2.7|) use 
has been made of the property [Q&, Uh] = 0, and the facts that the vacuum |0) is the only 
physical state in the PGM (at zero temperature) and that Uh\0) = |0) for any h. Therefore 
the expectation value (O)^ in the limit #ym — > is expressed as the average over h of the 
expectation values (0)p t h in the PGM with the boundary condition B h : 

{0) p ^ J dh±- jT VgVcDcDB O exp (i jT £ PGM ) = Jdh (O) , h . (2.8) 

The integrations over (c, c, B) should be done using periodic boundary condition. 

" We consider for simplicity a system without quark fields. The quark fields do not decouple from the 
(g, c,c, B) system if we impose a non-periodic boundary condition on g. 
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The PGM is still not easy to analyze for a general gauge-fixing function G[A^ c, c, B\. 
Fortunately the matter becomes remarkably simple if we adopt what is called the OSp(4/2) 
symmetric gauge |L5| given by the gauge-fixing function Gos P - 

Gos P = \ 5 A {tr (Al + 2icc) } , (2.9) 

where A is the gauge parameter and 5a is the anti-BRST transformation: 

S A g = -gc, 5 A c=^{c,c}, 6 A c=-iB, 5 A B = 0, (E = i{c,c} - B) . (2.10) 

This is because the action of the PGM is written in a manifestly OSp(4/2) symmetric form 
by introducing the superspace (x^,^,^) with Grassmann-odd coordinates 9 and 9 |Tl| : 

S OSp = jjd 4 x 5 A 5 B {tr (Al + 2icc) } = jjd 4 xj d9d~9 tr (ri MN d M G j d N G) , (2.11) 

where the superfield G(x, 9, 9) is defined by 

G{x,9,9) = (l + 95B + 95 A + 995 A SB)g(x) 

= g + 9gc-9gc + 99g(iB + cc), (2.12) 

and the superspace metric t]mn (M, N = 0, 1, 2, 3, 9, 9) is given as 

T]^ = diag (1,-1, -1,-1), rjg e = -r) el = i, others = 0. (2.13) 



OSp(A/2) is the rotation in the superspace (x^ } 9 } 9) which leaves the metric rj invariant. 
The BRST and anti-BRST transformation, 5b and 5a, correspond to the translation of 9 and 
9, respectively. 



Then, the Parisi-Sourlas dimensional reduction mechanism |T2[ tells that the PGM of eq. 
(|2.11| ) in four dimensions is "equivalent" to the chiral model in two dimensions described by 
the action 

S chira i = -y Jd 2 xtr (Wgidtf) . (2.14) 

The equivalence holds also in the present case of the finite temperature system in the real- 
time formalism with twisted boundary conditions Bh (|2.3| ) (the proof of ref. |16| which does 
not rely on perturbation theory applies to the present case): two spatial coordinates (x2,X3) 
and two Grassmann-odd coordinates (9, 9) in the 3+1 dimensional PGM cancel to leave the 
1+1 dimensional chiral model with the same boundary condition. The equivalence implies in 



particular that [11 



d 4 xe^%TN a (x)A b M)^ = -2n / ^^(^(^(O))^ 1 , (2.15) 
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where the LHS (RHS) is the Green's function in the 3+1 dimensional PGM (1+1 dimen- 
sional chiral model) with common /3 and the boundary condition Bh- In eq. Q2.15Q the four- 
momentum and the indices /i and v on the LHS should have components only in the two 
dimensional part /i, v = 0, 1 of the RHS. The reason why on the LHS of eq. Q2.15Q is 
converted to A^ = A^ — (A^ = gd^g) on the RHS is the relation = —iSa^b (AO- Note 
that A^ is the Noether current of the SU(N)r symmetry in the chiral model. 



From eqs. ( ]2.8|) and ( |2.15|) , we see that the confinement condition of KO fails if the chiral 
SU(N)r current A^ in the 1+1 dimensional chiral model contains a massless mode for h in a 
finite range of SU(N). At zero temperature {(3 = oo), the chiral model in 1+1 dimensions is 
realized in the disordered phase with a mass gap (the system is insensitive to the boundary 



condition Bh when (3 = oo) and hence the KO confinement condition holds [II]. We would like 
therefore to know whether the chiral model with a given boundary condition Bh undergoes 
a phase transition to an ordered phase having a massless Nambu-Goldstone mode coupled to 
A, 



3 Analysis of the 577(2) model 

Since we do not have a systematic non-perturbative technique to analyze the SU(N) chiral 
model in the real-time formalism with various boundary conditions Bh, we shall carry out 
the following approximate analysis to the SU(2) chiral model. Note that the SU{2) chiral 
model is equivalent to the 0(4) non-linear a-model via the expression g = cp l + iJ2l=i <£a& a 
((Pq + ip 2 a = 1). Rescaling ip^ and introducing the multiplier field a(x), let us consider the 
following 0(4) non- linear cr-model system: 

£o(4) = \&H> ■ - \a (V - 1) , (3.1) 

where the 4 component field (p = (<pi) = (<p , . . . , y? 3 ) is free from the constraint, and A is 
the bare coupling constant. The boundary condition Bh ( p.3|) for the 577(2) element h = 
exp (r] a a a /i) reads in terms of ip as 

B h : (pi(-T-ip,x) = (e T * Ta )..cp j (-T,x), (3.2) 

where the 4x4 matrix T a (a = 1, 2, 3) satisfies the property 

[T a ,T b ] = 2e abc T c , {T a 1 T b } = -25 ab l, (3.3) 
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and therefore we have 



e^^cosTyl + sinr^, (T,) 2 = -1, (r, = (r/ 2 ) 1/2 , T r , = Va T a /r/j . (3.4) 
The range of rj we should consider is < rj < 2n. 



Our analysis to the model (|3.1|) is to carry out the large-iV like calculation of the ef- 
fective potential to determine the phase of the model. Namely, we consider the effective 
potential which is the sum of the tree part and the one-loop (trace-log) term coming from 
the (^-integration. Since the large N expansion is a valid non-perturbative technique for the 
(ordinary) O(N) non- linear cx-model, it is expected that our analysis here will give a qualita- 
tively correct result for the present 0(4) model with a generalized boundary condition.[] Our 
question is whether the system ( |3.1|) with the boundary condition (|3.2|) undergoes a phase 
transition at f3 — f3 c {ff) from the disordered symmetric phase in the large f3 region to an or- 
dered phase with massless Nambu-Goldstone modes coupled to the Noether current (recall 
eq. ( |2.15| )). If at some f3 the 0(4) model is in the Nambu-Goldstone phase for a finite interval 
of the boundary condition parameter r) a , it implies that the KO confinement condition breaks 
down in the original vanishing ^ym limit of the Yang-Mills theory. 

It is straightforward to apply the path-integral formalism of ref. pi to the 0(4) model of 
eq. (|3.1|) with unusual boundary conditions Bh (|3.2| ). We need the Green's function Dp(x — y)ij 
(with the 0(4) indices i,j) on the contour of fig. [TJ satisfying 

(-□c - m 2 ) Dp(x - y)ij, = 5 c (x - y)5 ij , (3.5) 

and the boundary condition 

D^r-i^x)^ = (f« Ta ) ik Dp{r,x) kj . (3.6) 

In eq. ( |3.5| ) we have dc = (d/dr) 2 — (d/dx) 2 and 5c is the contour 5-function Q. Dp is 



obtained in the form 

D p {x -y) = D>{x- y) 9 c (r x - r y ) + D<(x - y) 6 c (r y - r x ), (3.7) 
with Dg given (in Fourier-transformed form with respect to the spatial variable) by 



£>!(t;A;i) = ;. - { ( 1 — e~ ptJ cos 77 ) 1 =f e~ p " sin 77 T, 



+ 



kM 1 



e Tlu}T 



cost] -e~ pu ')l± sin rjTr, e -P"^\ (3.8 



f3ui±iu)T 



" We do not know whether we can regard the present analysis as the N = 4 case of the large N expansion 



of an O(N) model since the boundary condition (3^) is particular to N — 4 and it breaks explicitly the 0(3)l 
symmetry. 
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where u = \ki\ and 



(3.9) 



For any boundary condition Bh, the generating functional Z[j] defined by the contour C of 
% 0, 

Z[j]= f V<p I Vaexpli f (C 0{ a) + 3 -<p + Ja)\, (3.10) 

JB h ./periodic I 7C V ' J 

factorizes in the limit T — > oo to the C1C2 and C3C4 parts, 



(3.11) 



We are interested in Z[j; C1C2] with infinite real-time segments. It is expressed as 
Z\j;C x C 2 ] = J J[ V<p A Va A exp{iJ™JtJdx(±v A -(Df) AB cp B 



where we have defined 



(3.12) 



3i(t, x) = j(t, x), j 2 (t, x) — —j It — 2^ ,X J ' 
and similarly for .1 A . Note that {ip t , (p 2 ) corresponds to (<p, cp) in TFD. 

■\AB 



(3.13) 



The finite-temperature propagator Dg (A, B = 1, 2) is given in momentum space as 



c s_ 

S+ c 



k 2 + ie 



( 1 
v 



k 2 + ie 




\ 

-1 
k 2 — ie / 



V 



2m8{k 2 )S + S. 
2m5(k 2 )CS + 



C S_ 

s + c 



-2m5(k 2 )CS_ 



2m8(k 2 )S + S- 



J 



(3.14) 



k 2 — ie 

where the 4x4 matrices C and S± are expressed in the form al + bT v and should satisfy 



C 2 - S + S_ 
S+S- 

cs± 



-P\ko\ 



{{cost] - e"^ 1 ^ 1 ) 1 - e(k ) sin r]^} , (3.15) 
.{(l-e^ fco cosr ? )l T e^ A;o sinr / Tje ± ^ fc o e ( A;o ) ) 

/ll fc o|) 



/(N) 

e -/3|fco| 
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with e(x) = sign (a;). We present here only the explicit expression of C: 



C = — 
>/ 2 /(l*b|) 

~ /3x COS 7] + 



«(N)i 



1/2 



e(ko) sin rj 



u(\k \) 



T 

J. n 



(3.16) 



where u(x) = (1 — e~ px cos rj + \J f{x) \ ' . Note that, in the particular cases of rj = and 
rj = 7r, the propagator ( |3.14| ) reduces to the familiar propagator for bosons and fermions, 
respectively. The inverse propagator (-Dg 1 ) appearing in eq. ( ft. 12] ) is given by 



D 



-1 



AB 



(ho, h] 



C -S- \f k 2 + ie 
-3+ c ) \ 
k 2 + i(C 2 + S + S-)e 
-2iCS+e 







-k 2 



ie 



C 

-S 4 



c 



-2iCS-t 
-k 2 + i{C 2 + S + SJ)e 



(3.17) 



Although the boundary conditions breaks explicitly the 0(3) ^ symmetry except the cases 
7] = and 7r, the effect of the breaking in the inverse propagator appears only at the ie parts. 

As stated at the beginning of this section, we consider the effective potential Vp((p A , a A ) 
which is the sum of the tree part and the trace-log term coming from the (^-integration: 



Vp(<p A ,a A ) 



1 



r 2 y* 2 ~ y 



1 

— 1— 
2 



d 2 k 
(2tt) 2 



tr In 



Vi 1 



ail 
-a 2 l 



(3.18) 



The phase of the theory is determined by the stationary condition of Vp((p A , a A ) with respect 
to ip A and a A (A = 1,2). The stationary points exist in the subspace act = a 2 (Vp (|3.18|) is 
singular when ot\ 7^ a 2 @). Therefore, the conditions to determine the groundstate are 







dcp j 



<x<Pa = °' 



0£l= CK2— Oi 



da a 



1 






N r J 


~ 2 

a.i =ao =a 




a ; 





dk 



-2PVk 2 +a 



AVk 2 + a f(y/W+ 



a 



(3.19) 
0, (3.20) 



where N = trl = 4 in the present case, and we have introduced the cut-off A for the k (spatial 
momentum) integration. The first condition fl3.19|) tells that we have either a = or ip A = 0, 
and from the second condition ( |3.20|) we see that <pl = ip 2 ,. The second condition (|3.20|) is 
rewritten using the mass gap m at zero temperature instead of Ao and A. Note that m is 
determined by eq. (|3.20|) with (3 = 00, a = m 2 and (p A = as 



2vr 



dk 



y/k 2 + m 2 



ln(^ 

V m 



(3.21) 
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which implies a familiar formula m = 2Aexp (— 2n/N\ ). Using eq. ( |3.21|) we have 



\A+1 9 v 

da a 



= \vl + = 0, (3.22) 

ot\=ct2=ct 



with K ^(a) defined by 



1 I _ p -2/3VW+^ 1 

h '^^l ^\-n^^n^=r--^T^\- (3 ' 23) 



Wk 2 + a f(^WTa~) Vk 



The property of our 0(4) model is qualitatively different between the periodic boundary 
condition (77 = 0) and the twisted one (0 < 77 < 2n). We shall discuss the two cases separately 
below. 

Periodic boundary condition 

It is easily seen that eq. ( ^.22| ) with rj = has no solution with a = for any (3 since 
Kp(a = 0) is positive and infrared divergent (note that f(x) = (1 — e _x ') 2 when 77 = 0). 
Therefore we have always <p A = 0. The expectation value of a, which is the (mass) 2 of the 
(p particles, is determined by eq. ( |3.22| ) with rj = and ip 2 = 0, and it is a monotonically 
increasing function of the temperature 1/(3. Therefore, the system with periodic boundary 
condition is in the disordered phase for any (3. 

Twisted boundary condition 

In this case, there is a critical temperature (3 c {rj) determined by the condition K^ = p o ^(0) = 
or explicitly, 

( m(3 c (rj) \ _ [°° d ( cosh(x/2)sinh(x/2) \ 

V 2 )~Jo dx{vnh 2 (x/2) + sm 2 (r,/2))- {6 ] 

The critical temperature (3 c {v) as a function of the boundary condition parameter rj (0 < 
r] < 2n) is depicted in fig. || in units of /3 c (7r) = 7re~ 7 /m (/3 C for the anti-periodic boundary 
condition (17, |ll|). As seen from fig. |2|, f3 c {v) is a monotonically increasing function of sin(r//2) 



and vanishes only when rj = 0. The meaning of f3 c {v) is th & t -^s(O) > (< 0) when (3 > (3 C 
((3 < (3 C ) (see fig. ||). The two phases separated by {3 c (v) are as follows: 



i) Low temperature region (3 > f3 c (7]): Eq. (|3.22|) has no solution of the type (cp A 7^ 0, a = 0) 



since we have Kp(a = 0) > when (3 > f3 c (see fig. |). Therefore, the system is in the 
disordered phase. The non- vanishing value of a determined by eq. ( p. 22 ) with ip A = 0, i.e., 



the intercept of the curve of fig. ^ with the a-axis, gives the (mass) 2 of the if excitation. 
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P C (T1)/P C (1C) 




. 2 



. 6 



. 8 



sin(r|/2) 



Figure 2: The critical temperature /3 c (r]) in units of (3 c {ir). 



ii) High temperature region (3 < /3 c (tj): Eq. ( |3.22j ) has no solution of the type ((f A = 0, a > 0) 
since Kp(oe), which is a monotonically decreasing function of a, is negative definite when a > 
(see fig. Therefore the system is realized in the ordered phase with (ip A ^ 0, a = 0). The 
Noether current of SU{2) R , A^ = (A») A=1 - {A^) A=2 (A» = e abc <f b d^ c + (p d^ a ), is coupled 
to the Nambu-Goldstone mode irrespectively of the direction of the expectation value ip A . 

The above results are summarized as follows. When (3 > /3 c (ir), our 0(4) model is in 
the disordered phase for all the boundary conditions Bh- When (3 < (3 c {tt), the models with 
t] in the range sin (f3~ 1 (f3)) < sin(j]/2) < 1 {f3~ x is the inverse function of (3 C ) are in the 
Nambu-Goldstone phase, while the models in the other range of boundary conditions are 
in the disordered phase. The range of the boundary conditions corresponding to the ordered 
phase increases as we increase the temperature 1/(3. Translating this result back to the original 
SU(2) PGM in 3+1 dimensions obtained as the vanishing gyyi limit of Yang-Mills theory, the 
KO confinement condition is satisfied at low temperature (3 > (3 c (tt), however it breaks down 
at high temperature (3 < (3 c {n) due to the contribution of the twisted boundary condition 
sectors in the average ( |2.8| ) with O = N ^A u . The Green's function ( |1.14| ) develops a massless 
pole when (3 < (3 c {n). 

Some comments are in order. First, we should comment on the compatibility of our result 



** The fc-integration in eq. (3.23) can be continued to the a < region, and eq. ( |3.22j ) with ip A = has 
a negative a solution even when /3 < f3 c . However, we do not adopt this solution since it implies that the (p 
excitation is tachyonic. 
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JKp(a) 







a 



Figure 3: The function Kp(a) with (3 below and above /3 c (r/). 



(i.e., that the Nambu-Goldstone phase is realized when (3 < (3 c (rj) for the twisted sectors) 
with Coleman's theorem which forbids the Nambu-Goldstone bosons in two dimensions. The 
origin of the absence of the Nambu-Goldstone bosons is that the (ordinary) massless scalar 
propagator in 1+1 dimensions 

(3.25) 



(2tt) 2 k 2 + ie 

does not exist because of infrared divergence. In the present case of the 0(4) model with 
a twisted boundary condition, the massless propagator in momentum space (see eq. (|3.14| )) 

the following form: 

1 



takes in the infrared k^ ~ 



(k) 



V 



k 2 



iire(ko)6(k 2 ) cot f|) T v 



VK 1 + COt 



(|) T„) e(k )5(k- 



-i7r(l-cot(^)r fl )e(feoW 



-P[l^ +i ne(k )5(k 2 )cot^) I) 



) 

(3.26) 

where V denotes the principal part. Coleman's theorem is evaded since the Fourier transform 
of the RHS of eq. ( |3.26| ) does exist: 

d 2 k 



(2n) 2 ^(^).iWW 



2n i gik-x 



\e{x% \e{x Q )9{x 2 \ 



(3.27) 



where 9{x) is the step function, 8(x) = (e(x) + l)/2. An intuitive reason why the deconfining 
transition occurs in our model is that the infrared singularity (in the perturbative ordered 
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phase) which caused the disordered phase at zero temperature is weakened by the twisted 
boundary conditions. The effect of the boundary condition becomes stronger as we raise the 
temperature and hence triggers the transition to the ordered phase. 

Our second comment is on the imaginary-time formalism. In this paper we have employed 
the real-time formalism of finite temperature field theory since our interest is in the KO 
confinement condition, which needs continuous four-momentum and cannot be discussed in 
the imaginary-time formalism. Forgetting this fact for the moment, let us consider what 
happens if we adopt the imaginary-time formalism defined by the straight vertical time contour 
[0, — i(3\ (cf. fig. |1|) in the above analysis of the 0(4) non-linear cr-model. Then, for the exactly 
periodic sector with r\ = we get the same conclusion that (a)n > and ((p)^ = for 
all (3. For the twisted boundary condition sector, however, we get a completely different 
result from the real-time formalism: the disordered phase with (y?)^ = is realized for all 
(3. This is because in the imaginary-time formalism with a twisted boundary condition, cp 
has no zero-mode to develop an expectation value. In particular, expanding ip into modes, 
cp(r, x) = X^-oo cp n (a;)e l ( 2n+1 ) 7rT / /3 (0<r</5) for the anti-periodic boundary condition sector 
(rj = 7r), the large N analysis shows that no modes cp n can develop expectation values. The 
expectation value (a) a is a monotonically decreasing function of the temperature 1/(3 and 
becomes negative for (3 < /5 c (7r) . Note that a negative (a)^ is not a trouble in this case since 
the effective (mass gap) 2 is given by (a)* + (tt/P) 2 , which is seen to be always positive. The 
discrepancy between real-time and imaginary-time formalisms in the case of twisted boundary 
conditions may seem strange, but it is not a problem since the twisted sector is not an ordinary 
statistical mechanics system. This discrepancy, however, is disappointing for an attempt to 
analyze the original chiral model by the Monte Carlo simulation since it is possible only in 
the imaginary-time formalism. 



4 Summary and discussion 

In this paper, as a first step toward the understanding of the deconfining transition in Yang- 
Mills theory in the sense of the breakdown of the color confinement condition of Kugo and 
Ojima, we have studied the model obtained by taking the limit of vanishing gauge coupling 
constant. This model at zero temperature is known to satisfy the KO confinement condi- 
tion. Adopting a special gauge- fixing function, the system in 3+1 dimensions in the real-time 
formalism is reduced to a "sum" of chiral models in 1+1 dimensions with various boundary 
conditions concerning the time contour. As a qualitative approximation to the SU(2) chiral 
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model we analyzed the equivalent 0(4) non-linear cr-model using the large- N like analysis. 
We found that the 0(4) model with a twisted boundary condition undergoes a transition to 
Nambu-Goldstone phase because the infrared singularity is softened by the boundary condi- 
tion. This implies the breakdown of the KO confinement condition in the 3+1 dimensional 
model. 

The pure-gauge model we considered in this paper, namely the vanishing #ym limit of the 
Yang-Mills theory, is physically trivial, and the transition we have found is not the singularity 
of the free energy for the 3+1 dimensional pure-gauge model. However, since the confinement 
mechanism of KO has an intimate relationship with the gauge field disorder in the direction of 



local gauge transformation ||14||, the breakdown of the KO confinement condition we observed 



should suggest the deconfining transition in the real Yang-Mills theory.^ 
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